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H ■ Abstract 
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Current constraints on photon velocity variability are summarized and displayed 
in terms of an energy dependent vacuum refraction index. It is shown that the 
energy-momentum balance of high energy Compton scattering is very sensitive to 

■ the outgoing photon speed. A missing energy observation in HERA Compton po- 
larimeter data indicates that photons with 12.7 GeV energy are moving faster than 
light by 5.1(1.4)mm/s. An asymmetry spectrum measured by the SLC longitudi- 



o 



nal polarimeter implies however an effect which is 42 times smaller, although the 
Cn . interpretation of the data is less clear here. 
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According to relativistic kinematics a photon velocity in vacuum c 7 does not 
depend on its energy u, while a possible dependency is constrained by the 
current photon mass limit m 7 <10 -16 eV [1] as 1 — c 7 (a;)/c < 10~ 32 uj~ 2 eV 2 , 
where c is a massless particle vacuum speed. However, the laboratory or stellar 
vacuum always contains background fields (matter) and quantum interactions 
can slow down or speed up photon propagation. Tiny changes of the pho- 
ton velocity have been predicted [2], [3] for such non-trivial, polarized vacua 
modified by electromagnetic or gravitational fields, temperature or boundary 
conditions within the perturbative quantum electrodynamics which allows to 
derive inverse relative velocities (vacuum refraction indices n = c/c 7 ) mainly 
for low energy uj ^ m (mis the electron mass) photons [4]. Even in the absence 
of background fields vacuum quantum fluctuations can influence light propa- 
gation as pointed out for the gravitational vacuum by recent developments in 
quantum gravity theory [5a-c]. Changes of photon speed are expected to be 
significant at photon energies close to the Planck mass ~ 10 19 GeV decreasing 
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with lower energies. Hypothetical Lorentz symmetry deformations considered 
for explaining the observed ultrahigh energy cosmic rays above the GZK cutoff 
(and possibly neutrino oscillations) [6] may also introduce an energy depen- 
dent photon speed [7]. 



2 Experimental Limits 



Magnitudes of these predicted effects are small and though may exceed by 
many orders the constraints imposed by the photon mass, all experimental 
tests so far show that different energy photons in vacuum move at the same 
velocity (light vacuum speed c) within the constraints displayed on fig.l (use 
of vacuum refraction index n{uj) instead of photon velocity is convenient to 
distinguish between photon mass and vacuum properties). 
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Fig. 1. Experimental constraints on the vacuum refraction index. 



The most stringent limits are coming from the detection of highest energy 
proton and 7 cosmic particles as first noted in [8] , since in a dispersive vacuum 
they would quickly decay by vacuum Cherenkov radiation p — > pj (n > 1) and 
pair creation 7 — > e + e~ {n < 1). These processes are kinematically forbidden 
in case 
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for Cherenkov radiation and pair creation respectively with M, E the proton 
mass and energy. Excluded areas in fig.l correspond to a highest detected 
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proton energy of E = 10 20 eV [9] and to a cosmic photon spectrum up to 
Umax — 22 TeV [10]. Also shown is a limit inferred from the highest observed 
electron energy of 2 TeV [11]. Other areas are excluded by experiments utiliz- 
ing direct time of flight techniques sensitive to \n — 1| Atc/D, where At is 
a time difference between arrivals of simultaneously emitted photons with dif- 
ferent energy and D is a distance to the source. While laboratory experiments 
are limited by time resolutions of typically a few psec and distances of a few 
km (an early SLAC result [12] \n — 1| < 2 • 10~ 7 is shown on fig.l by a narrow 
white bar at 15 GeV <uo <20 GeV) the astrophysical observations could do 
much better owing to huge distances to the source. In ref. [13] one can find 
limits on light speed variations in wide energy ranges based on different as- 
trophysical events; these limits suffer, however, from very uncertain distance 
scales. Meanwhile an observed spectacular gamma ray burst GRB990123 [14] 
followed by an optical counterpart detected within At = 22sec, with a distance 
z = 1.6, could establish a constraint \n — 1|<3- 10 -18 for 2 eV<uj<5 MeV, 
which is anyhow the order of constraints quoted in ref. [13]. Photons with high- 
est observed energies 0.35 TeV <uo < 10 TeV from a well defined active galaxy 
source (Markarian 421) put constraints |n — 1| < 2.5 ■ 10~ 17 c<j [15] (hatched 
area in fig.l). 



3 Compton scattering in dispersive vacuum 

Apart from the discussed threshold effects for vacuum Cherenkov and pair 
creation, the dispersive vacuum will modify the kinematics of other processes 
involving free photons according to the dispersion relation k 2 = u 2 (l — n 2 ). 
However, the tiny refraction imposed by such vacuum becomes observable 
only at high energies with corresponding small angles. When the photon (four- 
momentum k) interacts with a particle (four-momentum P) the vacuum index 
will contribute to the convolution Pk as 



where £,7 ^>1 are energy, Lorentz-factor of the particle, and 9 <C 1 is the 
angle between the photon and the particle. Thus, such processes in general 
could detect a relative photon speed variation, at given energy cu, as small as 
the order of I/27 2 . 

Below we concentrate on photon scattering off an ultrarelativistic electron and 
apply (2) in energy-momentum conservation to get sensitivity of the high en- 
ergy Compton process to the vacuum refraction index. If oj q , 9 , oj, 6 designate 
energy and angle of the incident and scattered photons, for ujq <^ u ^> m we 
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where j,£ are the Lorentz-factor and energy of the initial electron, 

4 7 cu sin 2 (gp/2) 
m 



, and n is the index for the direction 9 and energy c<j. In a case of laser Compton 
scattering on accelerator electrons the initial states (x, 7) are known to high 
degree of precision (typically to 0.01%) which allows to gain information about 
n from each event measuring the u and 9 (or the energy and angle of the 
scattered electron £', 9', since uj = £ — £', 9 = 9 '£' ' /u). Alternatively one could 
detect only the Compton edge i.e. maximal (minimal) energy of the scattered 
photons(electrons) u m = u at # = (£' m = £ — u m ) to measure n(u m ) down to 
values of 
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which follows from (3) if u m is measured with relative uncertainty Au m /u m . 
A dotted line in fig. 1 shows the potential of laboratory Compton scattering 
in limiting n according to (4) for optical lasers (uq ~ 2 eV) with a modest 
precision Au m /u m = 1% up to a photon energy of 100 GeV. 

The laser scattering is particularly attractive to test vacuum birefringence 
since the highest energy scattered photons preserve the laser polarization [16] 
which is easy to change. Flipping the laser linear polarization one could mea- 
sure ra_i_, ray components for multi-GeV photons by detecting the Compton 
edge dependence on _L, || polarization states (current bounds on the vacuum 
birefringence [17] are set by polarimetry of (near)optical photons coming from 
distant astronomical sources). 

In ref. [18] it has been proposed to test different quantities related to photon 
velocity by high energy Compton process measuring simultaneously the scat- 
tered photon and electron energies. However, this set of measurements is not 
sensitive to the photon speed which is accessible only from the photon energy 
and momentum combined information. To measure the photon momentum 
one has to register scattering angle of the photon or electron relying for the 
latter case on energy-momentum conservation. It is possible to indirectly reg- 
ister zero scattering angle of the photon by detecting the Compton edge as it 
pointed out above and only then an energy measurement alone is sufficient to 
obtain information about the photon speed. Another distinguished kinematic 
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point in the Compton process, where circularly polarized photons interact with 
longitudinally polarized electrons, is the energy asymmetry (between spin 1/2 
and 3/2 states) zero crossing which occurs at the maximal scattering angle 
of the electron and therefore at a fixed photon momentum. Thus, the corre- 
sponding energy u>a=o of the scattered photon gives a measure of the photon 
speed 
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1 = 



2 7 2 




- 1 



(5) 



This is a reduced form of eqn. (3) where angle detection is replaced by an energy 
measurement at the expense of dealing with polarized beams and is useful 
because most of the laboratory Compton devices are working as polarimeters. 

Derived relations allow to extract the refraction index and associated photon 
speed from existing polarimetric data. 



4 HERA polarimeter spectra analysis 



Consider photon spectra (fig. 2) from ref. [19] measured by the HERA Comp- 
ton polarimeter. The spectra were obtained by directing a CW 514.5 nm laser 
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Fig. 2. HERA polarimeter Compton events on top of background Bremsstrahlung 
and background subtracted Compton spectrum (inset) with fit results. Upper 
scale: original energy calibration using nominal Compton edge (GeV replaced by 
arb. units). Lower scale: recalibration using Bremsstrahlung edge. 

light against the HERA transversely polarized, 26.5 GeV electron beam with 
a vertical crossing angle of 3.1 mrad and detecting produced high energy 7- 
quanta with a sampling calorimeter. The whole detection scheme is designed 
for measurement of an up-down spatial asymmetry of the 7-quanta which is 
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introduced by a flip of the laser light helicity and is proportional to the elec- 
tron beam polarization while the energy measurement is auxiliary and serves 
as a mean to enhance the asymmetry by proper energy cuts. We are going 
to extract Compton 7's maximal energy from the spectra and estimate the 
refraction index via eqn.(3) at # = 0. Hence, following [19] and [20], we concen- 
trate on details of the experimental setup important for energy measurement 
only, ignoring all features related to polarization. 

The scattered Compton photons originate from an interaction region (IR) 
about 50 cm long, defined by the crossing angle and size of the electron and 
laser beams. Bending magnets downstream of the IR separate the electron 
and 7 beams and the photons leave the vacuum pipe through a 0.5 mm thick 
aluminum window to pass 39 m of air before entering the calorimeter which 
is installed 65 m away from the IR. 

Collimators placed at a distance of 47 m from the IR, define an aperture of 
±0.37mrad the same as angular size of the calorimeter as seen from the IR. 
The aperture is 15 times larger than the largest (horizontal) angular spread 
of electrons at the IR and 40 times larger than the characteristic radiation 
angle I/7 so the acceptance inefficiency can be ignored. The collimators are 
followed by magnets to sweep out charged background. 

The calorimeter consists of 12 layers of 6.2 mm thick tungsten and 2.6 mm 
thick scintillator plates surrounded by 4 wavelength shifters attached to 4 pho- 
tomultipliers. PMT signals from single photons are integrated within 100 ns 
gate then digitized with 12 bit ADCs and gains of the PMTs are adjusted to 
about 15 MeV per ADC channel. A fast DAQ handles the signals and operates 
without dead time up to an average data rate of 100 kHz. The detector perfor- 
mance has been simulated with EGS4 Monte-Carlo program and tested using 
DESY and CERN test beams. Measured energy resolution of 24% GeV 1/2 , 
spatial non-uniformity of ±1% and nonlinearity of 2% at 20 GeV are reported 
to be in agreement with the simulations. 

Apart from the laser light, the electron beam also interacts with residual 
gas, thermal photons and bending magnetic field in the beam pipe produc- 
ing respectively Bremsstrahlung, scattered blackbody and synchrotron radia- 
tion reaching the calorimeter. To measure this background the laser beam is 
blocked for 20 sec of each 1 min measurement cycle (light on/off is 40/20 sec). 
The procedure allows to eliminate the background by a simple subtraction of 
time normalized light-off spectrum from the light-on spectrum. Exact on/off 
durations are counted by DAQ clocks. 

At the time of the measurements an electron beam current of 0.32 mA and a 
laser power of 10 W provide 1 kHz rate above an energy threshold of 1.75 GeV 
while the background rate was 0.15 kHz. With such high threshold only the 
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Bremsstrahlung contributes to background since the scattered blackbody ra- 
diation maximal energy is 0.73 GeV and the synchrotron radiation is absorbed 
in the first tungsten plate of the calorimeter. 

Putting the laser photon, HERA electron energy and the crossing angle (u = 
2.41 eV, £ = 26.5 GeV, 0q = 7r + 3.1 mrad) into the definition of the kine- 
matic parameter x, we get x = 0.9783. The precision of the parameter is 
limited by the electron beam energy uncertainty a(S)/S ~ 10 -4 . Errors of the 
other constituents a(uj )/uJo ss 10~ 5 , a(m)/m 3 • 10~ 7 , A(9 ) ~2 mrad =>- 
A sin 2 (0 Q /2) w 3 • 10~ 6 contribute negligibly. 

To measure the ratio S/u m (the only unknown in the right part of (3) at = 0) 
we can utilize the Bremsstrahlung spectrum ([19], fig. 18) which helps to cancel 
the absolute energy calibration of the calorimeter since 

- = ^^ = ^ + 0(4. 10- 5 ) (6) 

uj m oiG m G m 



where a is a calibration constant and B m , C m are the Bremsstrahlung and 
Compton edges derived from the measured spectra in arbitrary units. It is 
easy to verify that influence of the term (2) to the Bremsstrahlung maximal 
energy is negligible i.e. a non-zero \n — 1| shifts only the Compton edge. 

A spectrum measured via calorimetry is conventionally described by a function 



where a parent energy distribution dT^/du incident on the detector is folded 
with a response function which is a gaussian with energy dependent width 
equal to the calorimeter energy resolution (in our case <t = 0.24 GeV 1 / 2 ), N is 
a normalization constant and E m is the cutoff energy of the parent distribution. 

The original energy calibration is made to match the nominal Compton edge 
Exj (1+x) = 13.10 GeV (fig. 2, upper scale) by applying a differentiation decon- 
volution method to find the cutoff energy. This method unfolds the spectrum 
by numerical differentiation to reveal a nearly gaussian peak (inverted) within 
the spectrum fall-off range and assigns the peak position to the cutoff value. 
The main drawback of this method comes from ignorance of the parent dis- 
tribution which results in a shifted answer in case of non-flat distributions as 
follows from (7). Therefore, to extract the B m , C m values from the spectra 
we have used a more precise approach (fitting via (7)) and have applied the 
differentiation method only to find the fit ranges around end points where the 
differentiated spectra peak, since outside of these ranges the spectra contain 
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no information about the cutoff energies. Such localization also helps to avoid 
possible bias of the fit results caused by physical effects affecting the spectra 
and not entering in the function F{E 1 ). Dominating among these effects are 
photon conversions between the interaction point and the calorimeter, detec- 
tor nonlinearity and spatial non-uniform response. These effects change the 
shape of the spectra in a way the function (7) is not able to describe ade- 
quately over the full energy range which is expressed also in ref. [19] and is 
noticeable for original fits shown on fig. 18 and fig. 21 of the same ref. 

Fitting the function F{E 1 ) to the background spectrum with the beam-gas 
Bremsstrahlung cross-section [21] as the parent distribution and 2 free, vari- 
able parameters E m , N, we get E m = 27.799 ± 0.047 = _B m , (see fig. 2). The fit 
range is predefined by numerical differentiation of the spectrum as discussed 
above. From a similar fit to the background subtracted Compton spectrum 
(from fig. 21 of ref. [19]) by F{E 1 ) with the Compton cross-section [22] as par- 
ent distribution, we find E m = 13.322 ± 0.010 = C m (fig. 2 inset). According to 
the derived numbers B m ,C m and relations (6), (3) we have the Compton edge at 
w m = 12.70±0.02 GeV, well below from the nominal u m {n = 1) = 13.10 GeV 
value and a vacuum index for the 12.7 GeV photons n — \ — (1.17±0.07)10~ n 
which is responsible for such reduction. 

Now we return to the above mentioned systematic effects to estimate their 
possible influence on the obtained cutoff energies. The non-evacuated path of 
7 beam line serves as an extended target to convert them into e + e~ pairs, sub- 
ject to continuous energy loss and multiple scattering before registration by the 
calorimeter. This modifies the spectra by enhancing lower energy parts with- 
out affecting the highest detected energies from non-converted 7-quanta. The 
most significant instrumental source affecting the result is the detector non- 
linear response Ea~ l (l + fE) under a given energy E with / = — 0.001 GeV' 1 
from the quoted nonlinearity of 2% at E = 20 GeV (/ < corresponds to 
a conventional calorimetric nonlinearity arising from shower leakage). This 
brings the ratio (6) to 

(l + f(£-u m )) (8) 



with a corresponding correction of 0.52 • 10 11 for n and half of that value as 
the correction error. 

Another possible source of the edges mismatch would arise if the Bremsstrahlung 
and Compton beams incident on calorimeter are separated in space. Propa- 
gating the quoted spatial non-uniformity of the calorimeter ±1% (B m /C m — > 
B m /C m {\ ± 0.01)) to the value of n we finally have 

n = 1 - (1.69 ± 0.07 ± 0.38 ± 0.26)10- n 
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with statistical and systematic non- uniformity, nonlinearity errors displayed 
separately. 

For completeness of the analysis, we discuss a few additional systematic sources 
which have no significant effect on the energy distributions. First is the ADC 
electronic pedestal with a width equal to ~ 30 MeV and a measured sys- 
tematic shift of the mean of ±4 MeV . This is neglected, since the pedestal 
spread is incorporated into the energy resolution while the shift is less than 
the result's smallest, statistical error by almost one order of magnitude. Next 
is an emission of multiple photons by a single electron bunch resulting in en- 
hanced maximal detected energies due to a pile-up in the calorimeter. Using a 
Poisson distribution and evaluating the quoted single photon emission prob- 
ability — 0.02 at the given Compton rate, one readily has probabilities 
for 2 photon emission pc(2) = 2.2 ■ 10~ 4 and ps(2) = 9.9 • 10~ 6 for Compton 
and Bremsstrahlung respectively. The latter number is too small to cause any 
considerable shift of B m , since the whole Bremsstrahlung spectrum contains 
less than a few pile-up events. Concerning the Compton edge, correcting for 
pile-up would only aggravate the observed energy reduction. The same is true 
also for non-linear Compton scattering events where an electron emits two or 
more photons at once. 



5 SLC polarimeter asymmetry analysis 

In ref. [23] one can find a Compton asymmetry (fig. 3) measured by the SLC 
polarimeter where high power laser pulses of 532 nm circular light interact 
with longitudinally polarized bunches of 45.6 GeV electrons under a crossing 
angle of 10 mrad and recoil electrons are registered by an array of Cherenkov 
counters installed downstream of two momentum analyzing bending magnets. 
Each channel of the detector integrates multiple electrons per pulse, within a 
certain energy range according to its position in the array. Following a detailed 
description of the polarimeter setup in [24] one infers that energies detected 
in the N-th channel are constrained by 

£Ln {max) = Co(S N + (-)D/2 + S- Sr 1 (9) 

where C = 296.45 GeV ■ cm, S = 10.58cm, Sn = N cm, D = 1 cm which is 
the channel size and S c is the Compton kinematic endpoint distance from the 
channel 7 inner edge which also depends on the initial electron beam position 
relative to the detector. Information about the photon speed is encoded into 
the relation of the Compton maximal and asymmetry zero crossing energies 
according to eqn.(5). A coarse granularity of the detector (binning in fig.3), 
however, makes it difficult to apply this simple kinematic method. Instead 
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Fig. 3. SLC polarimeter asymmetry (lower scale) with fit results (upper scale). The 
dotted line shows the parent distribution dY,\/dY, c . The lower part displays the fit 
residuals (right scale). 



one can utilize dynamic features of the Compton scattering in the case of 
n^l. Using an invariant representation of the Compton process in ref. [22], 
for longitudinal polarization of the incident electron beam one can write the 
cross-section as 

dE c dT^x nr 2 ( 1 , t . . \ 

+ A— — = — h 1 -y-4r(l - r) + Am (10) 



dy dy x 



where r e is the classical electron radius, A is the electron beam and circular 
light polarizations product, r — y/ (x—xy), u = rx(l — 2r)(2—y), y— 1 — Pk/Pko 
with ko being the photon's initial four-momentum, and x = 2PA; /m 2 , which 
is the kinematic parameter defined above. 

To introduce a refraction index into the cross-section, we modify Pk entering 
in y according to (2) and scale cross-section (10) by a factor of 

/ 2 dn\-i 
( U +niU d^) 

which accounts for a change of the delta function 5(lu 2 — k 2 ) to 5(n 2 u) 2 — k 2 ) 
in the phase space of the outgoing photon. In addition we use (3) and energy 
conservation to eliminate 9 and express the cross-section in terms of £' . 

The asymmetry An measured in a given detector channel iV is a product of 
A and an analyzing power I\/I c (A N = A • I\/I c ) where 
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(11) 



min 



with £' being the scattered electron energy limited by the channel's energy 



It follows from (9)- (11) that in the case of n — 1, the parameters S c and A 
establish horizontal and vertical scales respectively (energy and asymmetry) 
in fig. 3. However, these variables alone are not sufficient for a satisfactory 
description of the asymmetry distribution as indicated by a least squares fit 
performed with only two free parameters S c and A. Ref. [24] also reports 
about interchannel inconsistencies which dictate the choice and use of only 
one channel (number 7) for the polarization measurement. 

To extract the photon speed we add one more free parameter ip = 2 r y 2 (n — 1 ) 
and use the polarized Compton cross-section modified by dispersion, assuming 
a constant refraction index across the entire energy range of the measured 
asymmetry. Now the \ 2 minimization converges with A = 0.628 ± 0.009, S c = 
0.970 ± 0.037 and ip = -(6.49 ± 0.08) 10" 3 (fig.3), which yields 



for photons in the energy range 16.3 GeV < to < 28.3 GeV. 

An influence of the detector response on the asymmetry is quoted in [24] 
to be about 1%, which is much smaller than statistical fluctuations and we 
ignore it. Assuming perfect circular polarization of the laser light, A equals 
the electron beam polarization, which is measured to be 0.612 ± 0.014 from 
the channel 7 asymmetry. Both numbers agree within statistical and declared 
1.41% systematic [24] errors, and at the same time the n ^ 1 hypothesis 
allowed the asymmetry spectrum to be fitted successfully. 

Although the obtained result is more precise compared to the HERA obser- 
vation, it is less reliable because of the multi-electron generation-detection 
scheme and a theoretical uncertainty. The multi-particle mode, in general, 
poses difficulties to separate and treat the systematics and it also forced us 
to abandon the clear kinematic approach utilized in the case of the HERA 
polarimeter, while the method applied for modification of the Compton cross- 
section is somewhat heuristic and may introduce theoretical errors. 



n = 1 - (4.07 ±0.05)10 
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6 Discussion 



The observed value of the index, obtained from one sample of the HERA po- 
larimeter data, is statistically significant and does not contradict any previous 
experimental result (fig.l). It is below unity testifying that 12.7 GeV energy 
photons are moving faster than light (by c(l — n) = 5.07± 1.41 mm/s). How- 
ever a SLAC experiment shows that for photons of energy 16.3 — 28.3 GeV, 
the departure from the speed of light is at most 0.122 ± 0.0015 mm/s. 

Although the sign of the effect alone may be favorable for some theories dis- 
cussed in sec.l, the detected magnitude is too large to be associated with 
polarized electromagnetic or gravitational vacuum. So, the outcome is unex- 
pected, especially in view of the sharper limits for surrounding energies (see 
fig.l) and it is interesting to see whether the result can stand an examination 
by dedicated measurements and/or rigorous analysis of other pieces of data. 
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